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Abstract 

In this paper, we point ont and stndy a generic type of signals existing in the primordial 
nniverse models, which can be used to model-independently distinguish the inflation scenario 
from alternatives. These signals are generated by massive helds that function as standard 
clocks. The role of massive helds as standard clocks has been realized in previous works. 
Although the existence of such massive helds is generic, the previous realizations require 
sharp features to classically excite the oscillations of the massive clock helds. Here, we point 
out that the quantum huctuations of massive helds can actually serve the same purpose as 
the standard clocks. We show that they are also able to directly record the dehning property 
of the scenario type, namely, the scale factor of the primordial universe as a function of time 
a(t), but through shape-dependent oscillatory features in non-Gaussianities. Since quantum 
huctuating massive helds exist in any realistic primordial universe models, these quantum 
primordial standard clock signals are present in any inhation models, and should exist quite 
generally in alternative-to-inhation scenarios as well. However, the amplitude of such signals 
is very model-dependent. 
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1 Introduction 


In the past decades, through interactions between observations and theories, the Standard 
Model of cosmology has been established better than ever. With more abundant data coming 
from Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) in the near 
future, there are potentially different types of signals in data that could provide information 
beyond the Standard Model and allow us to address many deeper questions. In this paper 
we are interested in one of the most important questions in the early universe cosmology, 
namely, how to observationally distinguish the inflation scenario from the alternatives. 
Experimental data containing information on the primordial fluctuations determines how 
deeply we could address this question. Although many possible beyond-Standard-Model 
signals could provide invaluable information on model details within a scenario, unfortunately 
only very few of them could be used to model-independently distinguish different scenarios]^ 

A well-known example is the tensor mode of the primordial fluctuations [^, which may 
be observed as the primordial B-mode polarization in CMB |^. This signal records the 
magnitude of the Hubble rate in the primordial universe, and thus is able to at least model- 
independently distinguish inflation from some of the alternative scenarios with much more 
slowly varying scale factors]^ 

Another example that has been proposed recently is to make use of the oscillation of 
massive helds as the standard clocks to directly record the time-evolution of the scale 
factor of the primordial universe a{t) - the dehning property of a specihc scenario. In previous 
works, this idea has been realized in terms of the classical process [7|^. Namely, some kind of 
sharp features are necessary in the models to classically excite the oscillation of the massive 
clock helds. The ticks of the clock are recorded as special scale-dependent oscillatory features 
in the density perturbations, carrying direct information of a{t). Although the existence of 
the required massive helds should be general in any realistic model, the sharp features may 
not exist in any model despite of the fact that there are some well-motivated examples [^. 
We call this case the “classical primordial standard clock”. 

In this paper, we point out and explain that the same idea can actually be realized by 
quantum huctuations of these massive helds in the time-dependent background, which is 
now completely general. 

If the mass of a candidate clock held is large enough, there is always a regime during 
its evolution in which this mass is larger than both the physical wavenumber of the clock 
held and the mass scale of the horizon. We call this regime the classical regime. In the 
classical regime, the quantum-mechanically-originated oscillations of these massive helds 

^Nonetheless, we stress the current status that the inflation scenario is much better established than any 
other alternative scenarios. See for recent discussions. As we will see, the approach in our paper on this 
model-building aspect is phenomenological. 

^Having that said, some alternative scenarios still give similar prediction of gravitational waves. For ex¬ 
ample, both matter bounce and string gas cosmology predict nearly scale invariant spectrum of gravitational 
waves, and one needs to rely on details of the gravitational wave spectrum, as well as other observables, to 
distinguish such scenarios [^|^. 


2 



behave very much like those in the classical standard clock case, in the sense that these helds 
are approximately homogeneous over a region much larger than their Compton scales and 
they are oscillating as if in an approximately flat spacetime. For example, in the expansion 
scenarios, this classical regime happens in a relatively later stage during the evolution of 
the clock held, when the wavelength of the clock held gets stretched so that its physical 
wavenumber falls below its mass, and at the same time when the horizon scale is still below 
the mass; in the contraction scenarios, this classical regime happens in an earlier stage of 
the evolution, before both the wavelength of the clock held and the size of the horizon get 
contracted and hence their mass scales exceed the mass of the clock held. 

Once this is realized, the dynamics of recording the scale factor evolution a{t) becomes 
similar to that in the classical standard clock case. The ticks of these classical-like oscillations 
of the clock held get imprint in the correlation functions through the resonance mechanism & 
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However, unlike the classical standard clock case, the quantum huctuation of the massive 
held happens persistently and is not a source of scale-invariance breaking. Interestingly, this 
crucial diherence determines how the clock signals in the quantum case get subtly encoded in 
the density perturbations. They are encoded not as scale-dependent oscillatory features, but 
as shape-dependent osciWaioTj features in non-Gaussianities. We call this case the “quantum 
primordial standard clock”. 

For the inhation scenario, the ehects of quantum huctuations of massive held with mass 
of order the Hubble parameter H or larger on density perturbations have been studied to 
some extent 13 26 . The most striking feature is the distinctive scaling behavior of the non- 
Gaussianities as a function of momentum ratios. In the squeezed limit of the three-point 
correlation functions, the shape of non-Gaussianity S scales as [T^[T7l[T^, 


S oc 


’^long 


short 


Ps{cos9) , 


( 1 . 1 ) 


where 6 is the angle between the long and short mode (kiong and kghort) and s is the spin of 
the massive held. For this work, the most important feature in this formula is the scaling 
power a, which is determined by the mass m of the held. 
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The shape-dependent oscillation for the m > 3H/2 case is emphasized in 18 and interpreted 
as a consequence of the quantum interference between the inhaton and the massive held. 
In this paper, we point out another physical signihcance of this result - it is the quantum 
mechanical realization of the primordial standard clock and the signal directly encodes the 
scale factor evolution of the inhation scenario, a{t) ~ exp{Ht). This interpretation is impor¬ 
tant, because as mentioned, there are diherent kinds of beyond-Standard-Model signals from 
inhation models, and it is crucial to realize which ones can be used to model-independently 
distinguish the inhation scenario from the alternatives. We have briehy explained the rea- 
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sons behind this realization, and in the main text we shall illustrate them in details. For 
this purpose, it is also necessary to go beyond the inflation scenario and we shall study the 
effect of the quantum standard clocks in the alternative scenarios. We will show explicitly 
how their scale factor evolution a{t) is encoded in the clock signals. 

Massive fields are present in any UV-completed primordial universe models; they quan¬ 
tum fluctuate in any time-dependent backgrounds; and they couple to any other fields at least 
through gravity. Therefore, the signals of quantum primordial standard clocks are present 
in any inflation models, and should exist quite generally in the alternative-to-inflation sce¬ 
narios as well if the model-building would become more complete. However, the amplitude 
of such signals is highly model-dependent and its order of magnitude spans a wide range. 
In the gravitational coupling case, the amplitude is suppressed by the slow-roll parameters 
for inflation models [^; and if the mass of the clock held is much higher than the horizon 
mass, we expect additional Boltzmann suppression factors [^. On the other hand, for clock 
fields with mass of order the horizon scale (e.g. 0{H) in inflation) and for non-gravitational 
couplings, the signals may be potentially observable in the future [13 ,0)4 . The generic exis¬ 


tence of the signals and the variety of their amplitudes are very similar to the situations for 
the tensor mode. 

This paper is organized as follows. In Sec. we study quantum fluctuations of the mas¬ 
sive scalar held in arbitrary time-dependent background and emphasize the classical regime 
during which the massive field can serve as a standard clock. We use Fig. [^to qualitatively 
illustrate the relations of different scales in different scenarios, and the appearance of the 
classical regime in each scenario. In Sec. we first briefly review the mechanism of the 
classical standard clock, and then explain and summarize the main features and results of 
the quantum standard clock. We show how the scale factor evolution a{t) is directly encoded 
in the squeezed-limit bispectrum through an example. In Sec. and we show the details of 
the quantum standard clocks in the inflation scenario and the alternative scenarios, respec¬ 
tively. In Sec. we demonstrate through examples the highly model-dependent range of the 
amplitude of the clock signals in the inflation scenario. Sec. is conclusions and discussions. 


2 Quantum fluctuations of massive fleids 

For generality and simplicity, as in [^ we take the phenomenological approach and describe 
different kinds of time-dependent backgrounds of the primordial universe using a simple 
power-law function, 

a{t) = ao , (2.1) 

where oq and to are constants. This time-dependent background has a horizon size \t/{l—p)\. 
The requirement that the quantum fluctuations exit the horizon Axes the direction of t 
given the value of p. There are at least four interesting possibilities: |p| > 1 corresponds 
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to the fast-expansion scenario, namely inflation; 0 < p ~ 0{1) < 1 corresponds to the 
fast-contraction scenario; 0 < p 1 corresponds to the slow-contraction scenario; and 
— 1 <C p < 0 corresponds to the slow-expansion scenario]^ The time t rnns from 0 to -|-cxd 
for p > 1 and from —oo to 0 for all other p valnes. The conformal time r, which in this 
backgronnd is related to t by 

ar = , (2.2) 

1 — p 

always rnns from — cx) to 0. 

We now stndy the behavior of massive helds in this general backgronnd. We stndy scalar 
helds in this paper. Generalization to higher spin massive helds should also be interesting. 

The equation of motion of a massive held quantum huctuation, 6a, in arbitrary time- 
dependent background is 

6a + 3H6a — \di6a -I- m'^6a = 0 , (2.3) 

where dots denote derivatives with respect to time f; or in terms of the conformal time, 

6a'' + 2Ha6a' — dl6a + m^a^6a = 0 , (2.4) 

where primes denote derivatives with respect to the conformal time r. We will be interested 
in the cases where the mass is larger than the energy scale of the horizon, 

m > rriY, , (2.5) 


where 


mh = 


1 — p 

1 

t 

a r 


( 2 . 6 ) 


Even with the simple power-law (2.1), this diherential equation cannot be solved analyt¬ 


ically for arbitrary p. But the analytical expressions in the following limits are useful. 


2.1 Relativistic regime 

Denote the Fourier mode of 6a as Uk- Since the background is isotropic, Uk only depends on 
the magnitude of the mode k. When the physical wavelength of a mode is short enough to 


^We comment that, while the value of p is sufficient to define the scenario type, the terms “fast” and 
“slow” for different p-value ranges are named empirically. For example, in the available examples, the 
0 < p <C 1 scenario typically has very small Hubble parameter and the 0<p~D(l)<l scenario has 
large Hubble parameter, although this does not have to be the case just in terms of the phenomenological 


parameterization (2.1). 
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satisfy 


k/a^ m , 


(2.7) 


the EoM (2.4) is approximately 


v'l. + 2Hav'^ + k'^Vk = 0 . 


( 2 . 8 ) 


In the limit 


|fcr| Min(p, 1) 


(2.9) 


we can ignore higher order terms and the solution of this equation is approximately 


Vk —t 


-t .-1 / ^ \ ^-ikr 






( 2 . 10 ) 


where the normalization is hxed by the canonical quantization condition. We have also 
chosen the BD vacuum and we will further discuss this choice at the end of this section. 


Using (2.2) we can see that the condition (2.9) implies that the modes are subhorizon 


k/a 3> 1(1 —p)/t\. Nonetheless, for the massive held that satishes (2.5), the condition (2.7) 
is stronger than (|2.9). 


2.2 Classical regime 

When the physical wavelength of a mode is large enough to satisfy 

k/a m , (2.11) 


the EoM becomes approximately 

Vk + Sifhfc + rn?Vk = 0 . 


( 2 . 12 ) 


This equation can also be solved assuming (2.1), and we will be interested in the large |f| 
behavior. 


Vk \ — 


/ f \ -3p/2 

(-j (c+e-*™* + c.e*’"*) , 


(2.13) 


where the higher order terms can be neglected if 

\mt\ S> Max(p^, 1) . 


(2.14) 
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The parameter is the time at which k/a\t=t^ = m, and the quantization condition requires 

2ma^{tk) (|c+p — |c_p) = 1 . (2-15) 

To summarize, we define the classical regime to be the regime where the mass of the 
massive field is larger than both the physical wavenumber of the field and the horizon mass 


scale. The behavior (2.13) holds in the classical regime where both (2.11) and (2.14) are 


satisfied. The former condition makes the field oscillate homogeneously over a region much 
larger than the Compton scale ~ m~^. As we will see, m~^ is the typical time scale in the 
generation mechanism of the clock signals. The latter makes the massive field oscillate with 


frequency m as if in the approximately fiat spacetime. Note that the condition (2.5) is more 


general than (2.14) and is in principle enough for our purpose; in particular for the inflation 


case with p ^ 1, (2.14) is too restrictive and (2.5) is sufficient. Nonetheless, the condition 


(2.14) will be used for the alternative-to-infiation scenarios because it leads to a simpler 


approximation (2.13) 


This classical-like behavior reminds us of that of the classical standard clock, and is 
essential for the quantum fields working as standard clocks. Note that this is referred to as 
classical-like because the massive fields are not really created classically. 

In Fig. we sketch the evolution of several important scales in four scenarios. In the 
caption, we also emphasize several important properties that can be observed from this 
figure. 

As we can see from Fig. [T| for expansion scenarios, since the clock field originates from 
the relativistic {k/a > m > m^) regime, the BD vacuum is the natural initial condition 
to start these quantum fluctuations. The coefficients c+^_ are then determined by the BD 
vacuum. Nevertheless, as we will see later, although for the expansion scenarios we regard 
this BD state as the natural initial condition for the clock field, the most important property 
of the clock signal does not depend on this choice. 

For contraction scenarios, the clock field originates from the classical {k/a < m and 
m > TTih) regime. Only later on it evolves into the k/a > m regime which may or may not 
satisfy the condition m > mt (see Fig. [^. So the choice of (2.10) may not be natural at all, 
and the coefficients c+^_ are more arbitrary as long as the condition (2.15) is satisfiedj^ 


3 Quantum standard clocks 

3.1 Mechanism of the standard clock 

We start by briefly reviewing how the standard clock is able to record the time evolution 
of the background scale factor in the classical standard clock case. The classical os¬ 
cillatory behavior of massive fields is used as the standard clock. These massive fields, 

"^In fact, unlike the inflation scenario, since the contraction scenarios may or may not have an efficient 
way to get rid of the classical relics, even the classical origin of the clock field oscillation may be natural. 
Nonetheless, in this paper we focus on the minimal condition, namely the quantum fluctuations. 
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Figure 1: These figures illustrate the conditions for a massive field to be qualified as a quantum 
standard clock field in different scenarios. We qualitatively sketch the evolution of four scales (the 
constant mass of the massive field m, the physical wavenumber of the massive field /cdock/oj the 
mass of the horizon scale mhorizonj and the physical wavenumber of the massless curvature scalar 
mode fccurvature/a) as functions of the conformal time r or real time t in the inflation, fast (or 
slow) contraction, and slow expansion scenarios, respectively. Two different cases of m values are 
shown in some scenarios. The thick lines indicate the classical regime during which the massive 
field can be used as a standard clock. When the dashed lines intersect with the thick lines, the 
resonance between the curvature scalar mode and the clock field happens, and in all cases we find 
^ciock < /^curvature- At the resonance, fccurvature/a is always at the sub-horizon scales, fedock/a can 
be at either the sub- or super-horizon scales. 

















which are heavier than the mass scale of the horizon and homogeneous over a region much 
larger than 1/m, oscillate in standard ways in arbitrary time-dependent background. The 
most important property is its approximately constant frequency, ~ which can be 

regarded as a clock generating ticks. The ticks of the standard clock get imprinted in the 
density perturbations through the resonance mechanism |9^j^ This mechanism works in 
the following way. At subhorizon, the quantum fluctuations of the curvature scalar mode 
approach to those in the Minkowski limit. For example, if we choose the BD vacuum the 
mode oscillates as ~ in which K is the wavenumber. The physical frequency of the 

mode changes with time due to the expansion or contraction of the background. When this 
frequency coincides with that of the clock, resonance between the mode and the background 
happens, generating schematically the following contribution to the correlation function. 


^"^end 


^Tend 


"^begin 


dTg{t)e^'^^e g(t^) I dre 

^ ^Hoegin 


imt—iKr 


( 3 . 1 ) 


where g{t) represents some functions whose time dependence is much weaker than that in 
the other two factors, and [rbegin, Tend] denotes the period in which the resonance happens. 
The resonance happens if 


{mt — Kt) 
dr 


= 0 , 


(3.2) 


t=tf 


namely. 


a(f*) = K/m . 


( 3 . 3 ) 


The integration is proportional to 


P 

exp iimt^ — iKr^) = exp ( —i -m tiK/m) 

1— p 


( 3 . 4 ) 


where we have used the relation (2.2) for the power-law background. Notice that the function 


t{a) in (3.4) is the inverse function of a(f), so the scale factor evolution is directly recorded 


in the phase of this oscillatory signal and manifests as the relative spacing between the ticks 
of the signal in the iF-space. We call this signal the clock signal. 

The more precise formula that we often use is the following. 


^"^end 
'' ^begin 


1/2 //y\i/2p 


m 


\H, 


ko 


exp 


m 




,71 

l-pHk^\ko) ^^4 


p 


( 3 . 5 ) 


^In some special cases, the ticks may also be imprinted in a non-resonant way. We will discuss this in 


Sec. 4.1 and Appendix A.l 
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where ko is the mode that satisfies ko/a = m at some time to, and = p/to is the Hubble 
parameter at to- The =F applies when Hq is negative/positive. 

Away from the resonance point, mostly the integrand is oscillating rapidly with time and 
does not give significant contribution to the integration. At or outside the boundaries of 
this integral, especially at the Tend end, in some cases there may be a finite contribution; or 
in other cases, if the function g{t) is singular enough, the integrand can be divergent and 
needs to be cutoff or regulated by means of model-building which may leave some finite 
contribution to the integral too. In such cases, these contributions do not change the clock 
signal but add some overall shift. 

Previous examples of standard clocks make use of the classical clocks [^|^. They rely 
on some sharp process that excites the classical oscillation of a massive field. This process 
generates two different kinds of oscillatory features in the density perturbations, namely the 
sinusoidal sharp feature signal and resonance clock signal. These two kinds of features are 
smoothly connected to each other with specific patterns, and can be used as an identity 
for the standard clock signal in contrast to, for example, the resonance signals generated 
by the periodic features in the potential (See for a summary of other typical 

characteristics of the classical standard clocks.) The sinusoidal sharp feature signal is the 
one associated with the sharp feature. The resonance clock signal is the one induced by 
the subsequent oscillation of the massive field, and it records the time evolution of a(t) as 
illustrated above. 

There are at least two universal properties that make the clock signals “standard”: 

• Massive fields with mass larger than the horizon scale oscillate in standard ways in any 
time-dependent background. These include not only their most important oscillatory 
behavior, but also the typical time-dependence of the amplitudes which can also be 
used as auxiliary evidences. 

• The quantum fluctuations of the curvature scalar field, which are the source of the 
density perturbations, are at the subhorizon scales when the clock signal is generated. 
The behavior of such fluctuations approaches to that in the Minkowski limit and there¬ 
fore is universal for different primordial universe scenarios. This is very different from 
their highly model-dependent superhorizon behaviors. 


3.2 Quantum standard clocks 

In this paper we propose to use quantum fluctuations of the massive field as the stan¬ 
dard clock. In the time-dependent background, massive fields always fluctuate quantum- 
mechanically. We have noticed in Sec. |2.2 that, in the classical regime, the behavior of these 
quantum fluctuations resembles closely that of the classical oscillation. Besides sharing the 
above two universal properties, quantum standard clock exists in the absence of any sharp 
features. It only relies on the assumption of the existence of some heavy fields, which should 
be true in any realistic models. 
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To illustrate the mechanism and properties of the quantum standard clocks in arbitrary 
background, we use some simple examples of coupling terms between the curvature scalar 


perturbation ( (see 27 for dehnition) and the clock held perturbation 6a. For example, we 


use 


£3 ~ c^a^C^Sa (3.6) 

for the cubic coupling and 

C 2 ~ C2a?t6a (3.7) 


for the bilinear coupling, where 02,3 are some couplings whose time-dependence is weaker than 
1/m. The existence of these couplings require some non-minimal (i.e. non-gravitational) 
couplings between the inhaton and the massive held in the model; and such an inhation 
model example can be found in 13 . We discuss the gravitational coupling case in Sec. 

In most of this work, we do not pay much attention to the magnitude of the non- 
Gaussianity which we will later demonstrate to be very model-dependent. We concentrate 
on the momentum-shape dependence of the non-Gaussianity. 



Figure 2: An example of three-point function. The solid lines represent the curvature scalar mode, 
the dashed line represents the massive clock field. 


Using the in-in formalism (see 28,29 for reviews), these two coupling terms give the 


following contribution to the three-point correlation function of the curvature perturbation 
(see Fig.|^, 


(C')'D fdi, f dh{0\Hi{i,)ClHj{t,)\0y -2Re 

Jto Jto 

fO \ / rO 


tdh r dt2{0\CfHi{t,)Hi{t2)\0y 

Uto Jto 


= 2ul^Uk,Ukyr=o{ I df020^3^3 )( / drcsa^^u'yyy^ 

2,Zik^UkiUk2 |r =0 


(3.8) 


X 


dTiC2a'^VkoK^ / dT2C^a + / dTiCo,a / dT2C2a'^vl^u 


3/>,* 

kz 


(3.9) 


- 1 - c.c. -|- 2 perm. , 


where we denote the mode functions for C, and 6a as Uk and respectively; the prime 
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in (...)' means that we ignore the common factor (27r)^(5(^^ pj) in this expression. The 
“2 perm.” stands for two terms with the momentum permutation ki ^ and ^2 ^ ^ 3 , 
respectively. 


Let us hrst look at the term (3.8). The resonance between different mode functions only 


happens for the integral inside the second parenthesis. Analytical expression is possible if 
we look at the classical regime and the squeezed (^ 3 /^ 1,2 0) limit, in which the massive 

clock held has the wavenumber k^. In the classical regime, the massive held oscillates like 


the classical mode, described by (2.13). For the squeezed conhguration, the wavenumber 


of the massive mode is smaller than that of the two curvature scalar modes, and according 
to Fig. this oscillation provides a hxed frequency which at some point coincides with the 
frequency of the two shorter curvature scalar modes as their frequency decreases (in the 
expansion scenarios) or increases (in the contraction scenarios). This resonance happens 
when the curvature scalar modes Uk^ and Uk 2 are well within the horizon and described by 
their behavior in the Minkowski limit. 


Uk — 


—ikr 


i^/Aek 


(3.10) 


where e = We emphasize that the curvature scalar mode always starts in the 

subhorizon region quantum mechanically and then exits the horizon in all scenarios. The 
resonance always happens when the curvature scalar mode is still in the sub-horizon region. 


Therefore, the Minkowski behavior (3.10) is universal for all scenarios and this makes general 


analyses possible. On the other hand, although in this paper for dehniteness we always choose 
the BD vacuum for the curvature scalar mode, this choice is not essential for the clock signal 
and can be relaxed. Note that in the same squeezed limit the resonance does not happen 
for the other two terms with the momentum permutations. 

Note that the above conclusions only rely on the general properties in Sec. [^and Fig. 
and qualitatively summarize how the quantum clock signals are generated in general. 


Using these approximations, the integral in the second bracket of (3.8) becomes 


oc c_fci / dr I — 

“'''■begin 


-3p/2 


^ik\T 


(3.11) 


where we have ignored the coefficients that are independent of k^ and ki. The Tbegm and Tend 
denote the time period in which the massive held satishes the condition of being a quantum 
clock held, as illustrated in Fig. According to the convention in Sec. 2.2, is the time 
at which the curvature mode fcs satishes k^/a = m. 


Using the formula (3.5), we get 


c_ 


m 


\H, 


ki \ 


1/2 


2ki 

h 


1 

'2p 


exp 


m f 2ki\^^^ .71 

*1 -piLfcg V fca j ^*4 


(3.12) 
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In general, the coefficient c_ depends on and this introduces scenario- and model- 
dependent functions of in both the prefactor and the phase. The parameters Hk^ and k^ 
are evaluated at the reference point at which k^/a = m, and this reference point may be 
changed arbitrarily because remains a constant. For convenience we choose it to be 


hffcg and ^ 3 . The term (3.12) can then be written as 


m) (f 


1 

2p 


exp 


p 


m 


I-pH, 




— \ +i^{k^) 


(3.13) 


Note that / and p may also depend on p and m/Hk^. 

We can now see the reasons of our earlier comments that the BD state origin of either the 
curvature scalar held or the clock held is not essential for the clock signal. For the curvature 


scalar held, the non-BD component can also resonate with the clock held, similar to (3.11), 


although the component of the clock held that it resonates with is diherent. For the clock 
held, its non-BD state ahects the detailed fcs-dependence in c_, and is rehected as part of 


the arbitrariness of /(/cs) and 93 (^ 3 ) in (3.13), which is not the essential part of the clock 
signal. 


The integrals in (3.9) are more difficult to be carried out, and it will be one of the main 


goals of the next two sections where we work out these integrals in both the inhation and 
alternative-to-inhation scenarios. It turns out that, due to the classical-like nature of the 


clock held, the hnal result of the clock signal in the squeezed limit is proportional to (3.13). 


As a summary, the overall clock signal manifests as a shape-dependent oscillatory com¬ 
ponent in the squeezed limit of the three-point function. To standardize the convention, in 
this paper we will quote it as a component of the shape function S dehned in the 

following way ( 2 ^ 


{C^)^S{k,,k2,ks) 


ihhh 




(3.14) 


2=1 


where = 2.1 x 10 ® is the hducial power spectrum. In this convention, the clock signal 
in the squeezed limit kz/ki ^2 —)■ 0 is proportional to 


^clock ^ ^ 

h 

f2k. 


2^2p 


sm 


p- 


m 


2^ 2p 


sm 


u \ TW ) + ^(^3) 

p Pfc 3 V h 


m /2fciV^^ 

P - - TW + ¥^(^ 3 ) 


^h,fc3 \ ^3 


(3.16) 

(3.16) 


In the second line we have converted the Hubble parameter H to the horizon mass scale ruh 


dehned in (2.6). For non-exponential-inhation models these two parameters are diherent. 


and the horizon mass scale is more relevant physically. The parameter means that it 

is the value of ruh evaluated when k^/a = m. 
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We caution that, in quoting the above result from (3.8) and (3.9), we have assumed the 


prefactors Uki\T=o in those formulae to be proportional to l/y^. This is true for inflation 
and matter contraction; but may be non-trivial to realize for other scenarios. In any case 
this is only a convention. 


As we can see from (3.16), the most important clock signal is the oscillatory pattern 


sm 


P - - TW + ¥^(^ 3 ) 


^h,fc3 V ^3 


(3.17) 


Its envelop behavior also has some overall characteristic trend for each scenario, but may 
be subject to deformation due to model-building complexities]^ This clock signal may be 
extracted from data if we £x the long-mode and examine the bispectrum as a function 
of the short mode ki. This dependence gives an oscillatory signal that is similar to that of 
the classical standard clock . However, while the clock signal is a function of the overall 
scale in the classical standard clock case, here it is a function of the shape of bispectrum. 
These clock signals directly encode the evolution of a{t) ~ We also note that, with hxed 
fcs, for the expansion scenarios, the modes with smaller ki resonate hrst; for the contraction 
scenarios, the modes with larger ki resonate hrst. 

In Fig. 1^ for illustration and comparison, we sketch examples of clock signals (3.16) in 
four different scenarios. We have chosen the parameters m/mh,fc 3 so that the four examples 
have similar oscillatory resolution. The following is some comments on the values of m/mh^ks 
used in this hgure. Although for inhation ruh is constant, for alternative scenarios it is not 
and so becomes fc-dependent. Even if we choose the value of to be less than one for 

some scenarios in this hgure, it does not mean that the clock mass is less than the horizon 
mass when the clock signal is generated. For example, for the slow-contraction scenario, 
longer modes resonate later and the scale factor a evolves very slowly; so the condition 
k^/a = m is satished long after the clock signal is generated for the short mode ki. Since 
the horizon size shrinks fast in this scenario, m/ruh has dropped signihcantly to m/mh^fcg, 
even though it is larger than one when the clock signal for the short mode ki is generated. 
That is why in Fig. the value of m/mY^^ks is so small for this scenario. Vice versa for the 
slow-expansion scenario. 


® There is another important envelop factor that we ignore here, which is the analogous Boltzmann 
suppression factor that describes the probability of quantum mechanically creating a massive particle from 
the time-dependent background. For inflation this factor is of order ^ exp(—Trm/iJ) [^. Because H is 
approximately constant for inflation, this factor only gives an overall suppression to the amplitude, but not 
any fci/fca-dependent modulation. This is no longer true for the alternative-to-inflation scenarios. In these 
scenarios, we expect some analogous suppression factors too if the clock field oscillation has the quantum 
origin, with H replaced by rrih- However, since rrih is no longer a constant, these factors are expected to 
give some additional envelop behavior to the clock signal. These modulations are only negligible if we look 
at very short ranges of ki/k^. 

We also comment that for some of the alternative scenarios, classical origin of the clock field is also possible, 
and in this case the amplitude may not be exponentially suppressed. 

We leave these interesting issues for future investigations. 
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Figure 3: Examples of clock signals in four different scenarios: inflation (p = 20, m/m-h^feg = 5), 
fast contraction {p = 2/3, m/ruh^fca = 0.5), slow expansion {p = —0.2, m/ra\^^k^ = 2000), and slow 
contraction (p = 0.2, m/mh,fc 3 = 5 x 10“®). Also, see the last paragraph of Sec. 3.2 and footnote]^ 
for some comments on this figure. 


4 The inflation scenario 


In this section, we show the details of the three-point correlation function in Fig. in a 
single field inflation model with a clock field, which is essentially a quasi-single field inflation 
model. We pay special attention to the clock signal in this correlation function. 

Note that the purpose of this section is not meant to present a complete answer for 
a model. Instead, it only provides one example among several generic coupling terms in 
the quasi-single field inflation models (^, 14,Some of the results in this calculation 
have been obtained in 13 for m < 3H/2 and in [I^ 17 for m > 3H/2. Here we provide 


more complete and complementary results for this example. While most of the intermediate 
results apply to arbitrary m, we will be most interested in the mass range m > 3H/2 in the 
final results. 

For inflation, the equations of motion for the massless curvature scalar Uk and the massive 
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clock field Vk are, respectively, 


// ^ / I ; 2 n 

Uk - Uk + K 'Ufc = 0 , 

r 

„ 2 , ,9 

+ A: = 0 . 


(4.1) 

(4.2) 


Choosing the BD vacuum for both helds and normalizing the solutions using the canonical 
quantization conditions, we get 


H 


2^Mp A;3/2 


(1 + ikT)e 


—ikr 






(4.3) 


where 


fi = 


IP 


9 

4 


(4.4) 


In the definition of the shape function S (3.14), the formula for the power spectrum in this 
case is = P^/(87r^eM|). 

Since the single held inhation model is a well-established model, we are able to use 
the full mode functions (4.3) without having to use various approximations in Sec. and 
Sec. Some qualitative properties that we demonstrated in Sec. and Sec. such as 
the appearance of the classical regime and the functional form of the clock signals, will be 
borne out in the rigorous calculations. We will also see some additional properties that are 
not so obvious from the qualitative arguments, such as the appearance of the Boltzmann 
suppression factor and the clock signal from the non-resonance integral. 

We would like to compute two types of integrals - (3.8) which only involves single layer 


integrals without time-ordering, and (3.9) which includes two-layer time-ordered integrals. 


A summary is provided at the end of this section. 


4.1 The non-time-ordered integral 


Plugging in the mode functions (4.3), both integrals in (3.8) can be done exactly. The result 


IS 


+ c.c. + 2 perm. , 


(4.5) 
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where 
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T, — 

■y/F sec 


16^2 r(i 

/_ 

2(/ci -|- k2) 

1 

h 


—Ifl 


+i7r 


+ */i)r(l — i/i)r(—i/i) 2F1 ( - + 1 + 


2{ki + k2)y^ • ^ 17 ( 5 1 • fki + k2 

c.ch(.^)r(- - 1 - V, ( 


ki + /C2 


- 2 ' 


- 2 ' 


(4.6) 


In (4.5), “2 perm” stands for two terms with the momentum permutations ki <H- /C 3 and 


^2 -H- fcs, respectively. 

We can already see the shape-dependent oscillation factors [2(fci + k 2 )/as long as fi 
is positive. The clock signal is cleaner if we look at the limit in which the momentum triangle 
configuration is very squeezed, k^/ki <C 1 /y^, and hence the hypergeometric functions in 


(4.6) can be approximated as 1, 


very squeezed _ ^ , - 3/2 

-L\ 7 ^ fJjJ 


X 


. csch( 7 r/i)r(|-i/i) ^ 4 fciY^ 5 , . ^ f4:ki^ 

” - ■ - —'irj 


+0 


r(i — ip) 

— 2dii/i 


l-3/2 / h 

' Ua 


(4.7) 


Notice that, in the limit k^/ki <C 1/y^, the two momentum-permutation terms in (4.5) no 
longer contribute to the clock signal. 

Observationally we would be most interested in the region where the mass of the clock 
held is comparable to the Hubble parameter H. If the mass becomes much larger, the 
quantum creation of the massive held should be Boltzmann-suppressed [^. To see this 
Boltzmann suppression factor explicitly, we can further take the large mass limit p S> 1, 


s: 


clock very squeezed ^ C2C3 3/2 -27r;r f 


large mass 


->6 4 


32 


eHMl^ 


e ^ \ ^ 


-1/2 


4A;i 


— Ifl 


+ i 


4fci 

h 


4/r 


TT' 


3/2 


8 eHMl^ 


*"2^3 •gi -2-Kii (Y 


k-. 


-1/2 


sm 


, Aki TT 

+ 7 

h 4 


4 - c.c. 

(4.8) 

(4.9) 


As explained in Sec. this clock signal is a result of the resonance ehect between the 
BD vacuum behavior of the inhaton and the clock held oscillation e*™* in the classical 


regime. From the exact mode function (4.3), we can see that this classical regime appears 
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as the following late-time expansion, 




1-f coth(7r/i) I' •r(f/i) ^—kr 

r(ip -F 1) 


— i- 


—ifi 


. (4.10) 


Strictly speaking this leading term behavior is valid for kr —\/Ti- However, for kr —/i 
(i.e. k/a ^ m) the characteristic oscillations necessary for the clock held are already present. 


although the oscillation phase is not precisely described by (4.10) 


We note that such a resonance only contributes to the scaling power —iji of the momen¬ 


tum ratio, while in (4.8) we have i/r as well. This comes from a non-resonance contribution of 


the clock held contributed by the 1st term in (4.10), and turns out to be a special case of this 


single layer integral in the inhation model. We elaborate this in more details in Appendix 

EH 

4.2 The time-ordered integral 


To calculate (3.9), it is convenient to switch the inner and outer integrand of the hrst integral 


and rewrite (3.9) as 

/•O 


r ‘0 pT\ 


2 / dria / dT 2 a ^ ^ 'yfcg'Wfc* Wfc* / dT 2 a + c.c. -1- 2 perm. 

J —OO J T\ J —OO J —OO 

(4.11) 

Then both T 2 integrals can be calculated analytically. Inserting the result of the T 2 integrals, 
the three-point function now contains one integral over ri. 




C 2 C 3 

eHM^ 


{kik2k3Y^‘^l2 + c.c. -f 2 


perm. , 


(4.12) 


where 

- r 

= 2^Hankel + 2^2F2 = / dr (J^Hankel + *^2F2) 

J — OO 

*^Hankei and J^2F2 are dehned as 


(4.13) 


16cosh(7r/i) ^ 


(4.14) 
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-^ 2 F 2 = 


(-rfcg) {-rh) + (-rfca)) 2^2 

2fi — i 

2-3+*^-/fc7fc^r^r(i/i)e-’r^+*(^i+*^2)'r [—rk^) (“'^^ 3 ) + (“''"^ 3 )) 2-^2 


2fi + i 


(4.16) 


where 


113 

2 F 2 = 2^2 ( i/i + -, z/i + -; i/i + -, 2z/i + 1; 2zrfc3 
1 1 3 

2 F 2 = 2i^2 ( - - *)W, - - 2 “ 1 “ 2z/i; 2zr/c3 


(4.16) 


The J^Hankei term comes from the 2nd time-ordered integral in (4.11): the inner integral of 


this double integral has a constant part which is then turned to the -Tnankei term. 

In the remainder of this subsection, we shall show that Xnankei can be calculated exactly, 
and has clock signal in it. The X 2 F 2 can be calculated approximately overall, but exactly in 
some special limits; and it does not contain clock signal. 

• The Xnankei part. The Xnankei part can be integrated directly: 


-Hankel 


sech(7r/z) 

8^2 r(i/x + 1) {ki + k2fl‘^ 

x< T z/i ) [vr - zcosh(7r/i)r(z/xl)r(-z/i)] 2 X 1 


(4.17) 


where 


2 X 1 

2 X 1 


^ , 5 z/i 7 ill ^ 
2-^11 7 + y > ^ 1 


5 z/z 7 z/i 

2^1 1 7 - y>4 - y;i 


z/i; 


ki -h k2 


ki + k2 
k?. 


- 2 ' 


- 2 ' 


(4.18) 


Again, one can already observe the oscillatory factors [2(/c2 + ^ 2 ) 7 ^ 3 ]^*^ for positive /i. The 
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clock signal is cleaner once we expand (4.17) aronnd the very sqneezed limit fcs/fci \l 


very squeezed 


-Hankel 




71 r (ip + I) (csch(7rp) + sech(7rp)) /Aki 
r(ip + 1) V fcs 


—111 


2ir(|-ip)r(ip) 

g27r/i \ k^ 


+ 0 


L-3/2 ( h 

' Us 


— 2ii/i 


(4.19) 


Again note that, in the limit k^/ki <C l/^/p, the two momentnm-permntation terms in 
(4.12) no longer contribnte to the clock signal. 

• The X 2 F 2 part. We are unable to integrate the X 2 F 2 part directly. Nevertheless, one 


can use the method of resummation following Ref. 16 . Namely, we hrst Taylor-expand the 


hypergeometric functions in terms of r around r = 0; integrate each term in the Taylor- 
expansion with respect to r; and then try to resum the resulting series or approximate the 
sum by a hnite number of terms. 

Following this procedure, we hnd that, after integrating over r, the resulting series is 
in terms of powers of fcmassive/(fcmassiessi + fcmassiess 2 )- Tliis is also true for the other two 
momentum-permutation terms. So at least numerically this method can be used to compute 
most of the momentum conhgurations (except for the folded limit) to arbitrary precision, 
although in some cases the convergence is slow. 

To analytically investigate whether the X 2 F 2 part contains the clock signal, as explained 
we need to look at the case in which the massive held takes the long mode and two (of 
the three) massless modes are the short modes, ^massive "C fcmassiessi, 2 - To do that we only 
need to look at the integral X 2 F 2 in which the wavenumber of the massive held k^ is the 


long mode; the two other momentum-permutation terms in (4.12) are not needed because 


they correspond to the conhguration in which the massive mode is a short mode. The series 
converge fairly well in this limit, and we hnd that the hrst two terms in the series already 
give fairly good approximation to the full result. To reach the similar precision order in 


terms of k^/ki as in (4.7) and (4.19), we need take the squeezed limit of both terms. 


^ squeezed 

X 2 F 2 - > 




3k 


3/2 


4(4p2 + 1)^2 8(4p2 + Q)kl 


+ 0 


l - 3/2 ( kl 

U 3 


- 5 / 2 ' 


(4.20) 


where the hrst term comes from the 0th order Taylor expansion and the second term comes 
from the 1st order. Notice that, unlike the previous cases where the condition for the 
very squeezed limit is p-dependent, here the validity of (4.20) only requires k^/ki <C 1, 


independent of fi. See Fig. |^for some examples. This is because, as mentioned, the Taylor- 
series from the resummation method is a series in order of k^/ki, and in the several leading 
order terms that we pick up from the series, the hypergeometric functions take the forms 


such as 2 X 1 ( 2, |; 1 ± ip; ^ ) and 2 X 1 ( 2, |; 1 ± ip; ^ ). Unlike those in Xi and Xj 
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Figure 4: We demonstrate that both the squeezed-limit approximation (4.20) and the large ^ limit 


exact result (4.23) give good approximations to the full X 2 F 2 - The vertical axis is the corresponding 
contribution to the shape function, and the unit of Scontrib is 2c2C3/{eHMp). We used the isosceles 
configuration ki = k 2 - 


these hypergeometric functions have only one /x-dependent parameter and hence the validity 
condition for the small k^/ki expansion no longer relies on the value of /i. 

One can see clearly that there is no clock signal in X2F2- 

• Large fi limit of X2F2- In the large fi limit (but not necessarily squeezed limit), closed- 
form exact results can be obtained using the method of resummation . After one Taylor- 
expands the hypergeometric functions -7-2F2 and integrates each term with respect to r, new 
hypergoemetric functions arise of the form 2X1(2, 3 / 2 ; 1 ±x/i; k‘^/{ki -|- ^2)^) or 2X1(2, 5 / 2 ; 1 ± 
ifi] k‘l/{ki + ^2)^). Note that those hypergeometric functions can be expanded in their series 
form as 


2X1(0, 6; 1 ± ifi] z) 


{o)k{h)k 

^ {l±ifi)k k\ ’ 


( 4 . 21 ) 


where {a)k = r(a -|- k)/T{a) is the Pochhammer function. Note that the only term which 
does not have the 1 ± z/x factor in denominator is the fc = 0 term. Thus in the large fi limit 
only the k = term contributes. The result is simply 


2X1(0, &; 1 ± ifi] z) 


fl^OO 


> 1 . 


( 4 . 22 ) 


This greatly simplihes the calculation, and the series can then be resummed. We get 


fl^OC 

-‘-2F2 - > 


- AgA) 

4 (fci ^ 2 ) V 


-1- c.c. , 


( 4 . 23 ) 


where 


. ,113 2/^3 

fi = 3F211, - - ifi, --ifi]--ifi,i- 2ifi] j ’ 

. ,0 3 -3 .5 . , 2^3 

72 = 3X2 ( 2 , - - ifi, - - x/x; - - x/x, 2 - 2 x/x; - ^ 


( 4 . 24 ) 

(4.25) 
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Figure 5: A comparison between the squeezed limit approximation, (4.19) + (4.20), and the full 
numerical calculation of the bispectrum where the massive field is the long mode. The isosceles 
configuration ki = k 2 is used. The vertical axis is the corresponding contribution to the shape 
function, and the unit of Fcontrib is 2 c 2 Cs/{eHMp). The black solid line denotes the full analytical 
result, and the green dashed line denotes the Hankel part separately. The blue dots line represents 
the numerical calculation. One observes that for small //, the Hankel part gives dominate contri¬ 
bution. For larger /r, though the Hankel part still contribute all the clock signals, the 2 T 2 part 
contributes a significant shift of the amplitude. 
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/s = 4 F 3 ( 2, 2, - - i/r, - - i/i; 1, - - i/r, 2 - 2i/r; 


2fc, 


h + h 


(4.26) 


Note that the “c.c.” in (4.23) is in addition to that in (4.12) 


We also compare this large fi expression with the results with hnite fi in Fig. One can 
see that it is a good approximation even at the smaller fi values. 

• ^^Hankei + 1.2Y2- The following are a few aspects of the net result Xnankei + X 2 F 2 - 


Comparisons between the squeezed limit approximation, (4.19) + (4.20), and the full 


numerical calculation where the massive held is the long mode are shown in Fig. 


Parallel to what we did in (4.9), we can further take the large mass limit of these results. 


The clock signal is again given by (4.19) in these limits, 
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I-1 large mass 


TT' 


3/2 


32 


pF-f ^2^3 3/2 

eHMf 


C2C3 


__,,3/2 -TT^i 

16 eHMf 


-1/2 


Aki 

h 


— 111 


+ c.c. 



-1/2 

Aki 

dvr 

k) 

sin 

/iln — 

h 

+ 

1 _ 


(4.27) 

(4.28) 


where we can see the Boltzmann suppression factor. 

Note that in these limits, the contribution that does not contain the clock signal, namely 
(4.20), is not suppressed by the Boltzmann factor but only by the power law ~ /i“^. The 


physical interpretation 18 is that such a contribution does not describe the quantum creation 


of massive particles from the inhationary background, and is captured in the low energy 
ehective theory. However the clock signal corresponds to the quantum creation of massive 
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particles and hence has to be suppressed by a Boltzmann factor; and this contribution is not 
captured in the low energy effective theory language. 


4.3 Summary 


To summarize, in this section we studied the quantum standard clock signal in the inflation 
model. We have presented the analytical and numerical results on the bispectrum term 
Fig. [^for a single held inhation model with a clock held. We are most interested in the clock 
signals in which there are shape-dependent oscillations. 

For the clock signal, we are able to obtain the exact analytical result. In terms of the 


shape function S (3.14), this is given by 


^clock _ 


C2C3 

eHM^ 


(Xi + Xnankei) + c.c. 2 perm. 


(4.29) 


where Xi and Xnankei are given by ( |4.6| ) and (4.17), respectively. 

In the very squeezed limit, k^/ki the appearance of the clock signal becomes 


cleaner, and they are given by (4.7)-|-(4.19). Note that the other two momentum-permutation 


terms in (4.29) do not contribute to the clock signal in this limit. 


Further taking the large mass m S> 77 limit, we can see the Boltzmann suppression 


factor. The result is (4.9)-|-(4.28). Comparing (4.9) and (4.28), we can see that (4.28) is the 


leading behavior in this limit. 


As we can see explicitly in, for example, (4.9) and (4.28), these clock signals are the 


special case of the general clock signal (3.16) in the p —)■ oo limit. As expected, the function 


of the momentum ratio in the phase of this oscillatory signal is the inverse function of the 
inflationary scale factor a ~ 

There is also another contribution X 2 F 2 , which we show does not contain any clock signal. 
For this contribution, in general we can only obtain a series expansion but not a full closed- 
form result. Nonetheless, the analytical result in the squeezed approximation {k^/ki 1) 


in which the massive field takes the long mode is obtained in (4.20); and the closed-form 


analytical result in the large /i, but not necessarily squeezed, limit is obtained in (4.23) using 
the method of resummation. 

In Fig. we show several examples of the complete numerical results of the shape 
functions. These shape functions contain both the clock signals and non-clock signals. We 
choose fi to be of order one in this hgure. 


5 The alternative scenarios (0 < 7 ? < 1) 

In this section, we study the clock signals generated in the alternative-to-inflation scenarios 
with 0 < p < 1. This includes the fast-contraction scenario 0<p~C>(l)<l, such as 
the matter contraction scenario [^, and the slow-contraction scenario 0 < p 1, such as 
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Figure 6: The numerical results of the full shape functions with /i of order 0(1) in the inflation 
scenario. These numerical results include both the clock signals and the non-clock signals. All 
momentum permutations are included in these results. The method of Wick rotation [13] is applied 
in the numerical integrals to increase the convergence speed in the UV. To show more clearly the 
oscillation patterns of the clock signals, we have multiplied ki/k^ to the shape function. The 
vertical gray dashed lines denote the peaks of the oscillation. The periodicity in log(A;i/A: 3 ) is 1/n- 
This agrees with our analytical result. 


the Ekpyrosis scenario 31 . As the inflation example, we assume one massless held that is 


responsible for the curvature perturbation and one massive held as the clock held. 

However, comments are in order before we proceed to the more detailed analyses. 

As we know, alternative scenarios are not well established as inhation. For example, 
in alternative scenarios, the assumption of single massless held typically does not lead to 
viable models; even with one massless mode, its superhorizon behavior may be highly model- 
dependent. 

Fortunately, the standard clock signals are largely insensitive to such details. They 
manifest as some fast-oscillating signals added on top of the density perturbations with 
otherwise much smoother scale or shape dependence. Equally importantly, in the regime 
where the clock signals are generated, mode functions exhibit some universal behaviors. 
This makes the general analyses possible. Therefore, although unlike the inhation case 
where we have the complete mode functions, in this section we can still proceed by utilizing 
these universal behaviors. We start by briehy reviewing these behaviors in the following. 
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For arbitrary p, the scale factor and the relation between t and r are as follows, 




P 

1-p 


T = 


hti-p 


pao 


(5.1) 


As we showed, the clock signal is generated when two cnrvatnre scalars are the short modes 
and at the snb-horizon region, and the clock held is the long mode and enters the classical 
regime. We thus approximate the two short curvature scalars using the universal sub-horizon 
behavior. 


Uk —t 


—ikr 


(6.2) 


a\/Aek 

where we have again chosen the BD vacuum; and we approximate the clock held huctuation 


5(t by the solution in the classical regime (2.13), 


1 ^p 


/ro\ 2(i-p) 

c+e -(-c_e 


\T / 



(5.3) 


In the classical regime the physical wavenumber of the clock held has ehectively been taken to 
zero because it is much smaller than the typical frequency scale of the resonance phenomenon 
~ m, or equivalently as we noted previously, the held is homogeneous over a region much 


larger than the Compton scale. In (5.3), the hrst component is the positive frequency mode. 


the second is the negative frequency mode, Tq is a negative constant, and 

V = ml(-aoT„y' 


( 5 . 4 ) 


is the ratio of the clock held mass to the horizon scale mass at tq. Although the horizon mass 
is usually time-dependent, here v is evaluated at tq and by dehnition a constant. 


-ar 


This ratio is an approximate generalization of p (4.4) dehned in the inhation case. 
For simplicity, we choose Tq = — 1 and denote 


Uk ^ 


—ikr 


Vk OC 


(—r) i-p 
1 

(—t) - 




(5.5) 


The pre-factors of fcj’s may be easily restored at the end of the calculation. We will also use 
the above Uk for the long curvature scalar mode, because for contraction scenarios the long 
mode is still at sub horizon during the resonance in most cases, as one can infer from Fig. frh 


^There is also a small corner of parameter space in which the resonance happens when the long mode 
is at superhorizon. In such a case, oc (—t)“^. As one can check, as long as 2 {i-p) < < P^-p) ’ 


the conclusions in this section remain the same and we get the same result (5.13) with a different constant 
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5.1 The non-time-ordered integral 


The non-time-ordered integral ( |3.8 ) is analysed in Sec. 3.2 Only the following term contains 
resonance in the squeezed limit k^/ki —)■ 0, 


cl / dr(-r)“ . 


oc 


(6.6) 


Note that only the negative energy mode Vk contributes to the clock signal - unlike the 
inflation case, the integrals that do not resonate generally do not contribute to the clock 
signals, and we demonstrate this observation in Appendix A.2|P1 


5.2 The time-ordered integral 


Using (5.5), the integrals (3.9) become 


3p 


oc 


'0 


dxix^ -t- 
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r^xi 


r-oo 3p / 1 
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'0 " V 
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r^p 




X / dx2xf-^^ ( 
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' Q-ik3X2 


(5.7) 


(5.8) 


where we have redefined Xi = —Ti and X 2 = —r 2 , and in (5.7) switched the inner and outer 


integrand for convenience. The massive field is in the classical regime, and for the same 


reason we stated below (5.3) we can approximate fcs = 0. The inner integrals can then be 
done analytically. 


We first work on the integral (5.7). Integrating over X 2 , we get 

2-p 


(-1+P)X2 


2(l-p) 


c+Ee -i(l - p)ux 2 ^ + c_Ee i{l - p)vxl 
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> 1-P 
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UX^ ^>1 I 2 ( 13 ^ 
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i{l-p)uxl ^ 


1 

-P 


- ^r(l - |) [-C+ (-f(l - p)iyy^^ + c_ (i(l - p)iy) 


p/2 


(5.9) 


where En{z) is the exponential integral function. In the 2nd line, we have expanded the upper 


prefactor. _ 

®As we demonstrate in Appendix A.2 the conclusion, that the clock signals are only contained in the 
resonant integrals for the alternative scenarios with 0 < p < 1, are made empirically through examining 
several representative examples, and a general study is currently unavailable. 
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limit of integral in the large z/a;| ^ limit, because the resonance happens when m|ti| S> 1 


(see (2.14)), namely, when uxl ^ is large. The lower limit of the integral gives a constant 


term in the third line. 


Combining (5.9) with the integrand in the outer integral and using the observation from 


Appendix A.2 we can identify all the terms that contain resonance and hence potentially 
the clock signals. 


c+c_- 


dxix^ 


2p 


l-P 2i{l-p)ux^ ^ -2ikixi 


+ ^r(l-|) c+c*_ - p)u) 


p/2 


|c_|^ (i(l — p)^) 


p/2 


dxiXi 


3p — 

2(i-p) i(l-p)raj —2ikixi 


( 6 . 10 ) 


Note that the 2nd term is due to the constant term in (5.9). 


We now work on the second integral (5.8). Similarly, the inner integral can be done 

1 -r 1 

analytically and we then expand it in the large vxl~^ limit, 

2-p 


-1+P)X2 


2(l-p) 


c^Ez i{l - p)ux2 ^ + c*_Ee -i(l - p) 


jUX. 


1-p 


xi 
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-^--Xi 
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* -i(l-p)raj 
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(5.11) 


Note that, at the inhnity end, the inner integral gives zero because x^^^ decays as X 2 —t oo 
for 0 < p < l|^ Combining with the integrand in the outer integral, we identify all terms 
that contain resonance as 


2p 


c+c_- 


u 


dxix^ g-2it.r. 


(5.12) 


It is interesting to notice that the terms with doubled clock frequency in (5.10) and (5.12) 


cancel each other, and the hnal clock signal is contained in the following term. 


^r(i-|) c+ch (-i(i-p)z/) 


p/2 


|c_p (i(l — p)z/) 


p/2 


dxiXi 


3p 


2(l-p) i(l—p)i/x/ ^ —2ikixi 


(5.13) 


This integral is of the same form as (5.6). 


We note that both (5.6) and (5.13) are non-vanishing only in the presence of the negative- 
energy component of the clock held. This is due to our choice of the BD vacuum for the 


massless curvature scalar in (5.2). Although the BD vacuum is the natural choice for inhation 


®For the inflation scenario (|p| > 1) and the slow-expansion scenario (—1 p < 0), this is no longer true. 
These cases need to be treated separately. The inflation case is studied in Sec. 
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models, as discussed this does not have to be the case for contraction scenarios. In the 
presence of the non-BD component, the integral terms that contain resonance would be 
different and the positive-energy component of the clock held can also contribute. 


Both integrals, (5.6) and (5.13), are divergent at r —)■ 0 and needs to be cut off or regu¬ 


larized through model-building. However, as we emphasized the advantage of the standard 
clock is that such details do not affect the essential features of the clock signal contained in 


these integrals. As in (3.5), using the saddle point approximation, we can extract the clock 


signal from these integrals. 


oc ( 2 fci)" 2 +^ sin 


p u p ( 


-2kiToY^^ + phase 


(6.14) 


where we have restored some pre-factors of ki and the factor of tq. This result matches with 


the general result (3.16) with 0 < p < 1. To see the correspondence between the constants 


in these two formulae, we assume that the wavenumber k^ satishes k^/a = m at r = tq. 
We then have the following two relations, k^/ao = m, = a'/o?\rQ = p/ ((1 — p)ao'ro). 
Together with (5.4), these relations connect the conventions in the two formulae. Finally, 


as we cautioned below (3.16), we note again that, when quoting the convention S here. 


Mfci|r-s>o oc 1/y/^ has been assumed for those prefactors in (3.8) and (3.9), which may be 
nontrivial to achieve for alternative-to-inflation scenarios. 

It is interesting to observe that, in the time-ordered double integrals, naively there are 


terms in (5.10) and (5.12) containing clock signals with doubled oscillatory frequency; but 


they get cancelled in the end. A similar cancellation also appears in the inflation case, where 
we showed that the integral X 2 F 2 does not contain clock signals. In the inflation case, this 
fact can be readily understood (^: as the wavenumber of the massive held approaches zero, 
two vertices in Fig. are separated by the lightcone and so the non-analytic terms in the 
propagator between them do not depend on the time-ordering. Thus the overall integral that 
contains non-analytic terms in k factorizes. Since we used the exact mode functions in the 
inhation example, this property is automatically present through the /c-dependence of the 
coefficients c+ and c_. But in the analyses of the alternative-to-inhation scenarios, we did 
not use any properties of c+ and c_ nonetheless we still see a result that indicates some 


kind of factorizability in the integral. So, we expect a different and more general explanation: 
when we set ^3 —?• 0 and use the behavior of the massive held in the classical regime, the 
massive mode at the hrst vertex of Fig. acts much like the background oscillation in the 
case of the two-point function of the classical standard clock. 


5.3 Summary 

To summarize, in this section we studied the quantum standard clock signal in the alternative- 
to-inhation scenarios (0 < p < 1 ) including both the slow and fast contraction scenarios, 

^*^Even if the complete model and mode functions were given, in the classical regime the clock field is 
mostly at the subhorizon scale for the contraction scenarios (see Fig. [T|). 
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again using the example Fig. Despite of the uncertainties and immaturity of the actual 
model-building for these scenarios, we explained how we can use the universal behaviors of 


the mode functions (5.2) and (5.3) to extract the essential features of the clock signal from 


the integrals. We showed that both the non-time-ordered integral (3.8) and time-ordered 


integral (3.9) are proportional to the same resonant integral, which gives rise to the clock 


signal (5.14). As expected, the functional form of the phase of this oscillatory signal is the 
inverse function of the background scale factor evolution, a ~ t^. 


6 Amplitude of the signals 

There are a variety of well-established inflation models and it is possible to introduce some 
massive helds and study them case by case. In this section, we briefly discuss the amplitude 
of the clock signals using a concrete model of quasi-single held inhation 13 . We leave the 
alternative-to-inhation models to future study. 

In this model, besides the universal gravitational coupling, the inhaton is coupled to a 
massive held with mass of order H through a direct derivative coupling, and the massive 
held also has non-linear self-couplings. Following the notations in [^, we denote the sizes 
of these couplings as Oq/H and V"/H (for the cubic coupling), respectively. Both couplings 
are required to be much less than one due to the limitation of the perturbative calculation. 
The natural embedding of this model in the supergravity theory and its generalization are 
shown in [Tl]. The amplitude of the three-point function due to a variety of couplings have 


been estimated in the Appendix A of 13 and in (l4j. We comment that, in the diagram 


that contains the massive held self-coupling, the clock signal is still generated in the same 
qualitative way between the two short massless modes and one long massive mode, as we 
described in Sec. and The intermediate propagators of short massive modes do not spoil 
the resonance, because the time integration contains a region in which the phases of the 
massive mode and its complex conjugate cancel each other in the propagator. We can also 


directly use the amplitude estimates in 13,14 to get some ideas on the range of the size of 


the clock signals, since it is a smooth analytic continuation across the m = 3H/2 point. 

• If the direct coupling between the inhaton and the clock held is turned oh and there is 
only the universal gravitational coupling, the clock-held-mediated non-Gaussian cur¬ 
vature scalar correlations appear as loop diagrams, since the 6a^s always appear in 
pairs in the vertices of the gravitational coupling. The amplitude of these correlation 
functions are suppressed by factors of e and ~ 10“®. For example, using the cubic 


and quartic couplings in 32 , we can estimate that for the bispectrum /nl is propor¬ 


tional to eP^ with some possible enhancement factors from the UV cutoh scale, which 
may compensate the P^ suppression in non-Gaussianity. 

Gonsidering the direct coupling between the inhaton and the clock held through the 
derivative coupling turns on a bilinear coupling between them. Depending on whether 
the three-point vertex is of the type or (6a‘^ and the derivative structures in 
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these vertices, the amplitude of the bispectrum cau be /nl oc 0{e{9/HY), 0{e), or 
Here e is 0(0.01) or much less, aud {0/HY -C 1 for perturbative calcula- 

tiou. 


• Cousideriug self-iuteractious of the clock held turus ou a larger source of uou-liuearity, 
aud the resultiug bispectrum has au amplitude 

/nl oc p-^^\v'"/H){eo/Hf ~ l9p{y'"/H){eo/Hf , (6.1) 


agaiu V" / H aud 9q/H are much less thau oue. This possibility gives the largest uou- 
Gaussiauity amoug the examples meutioued here. Note that the amplitude of /nl here 
is evaluated at the equilateral couhguratiou. For the squeezed couhguratious, we also 
have to take iuto accouut the euvelop behavior iu (3.16) aud au extra factor of k^/ki 
due to the reductiou of couhguratiou space wheu detectability is coucerued, both are 
suppressiou factors for the iuhatiou sceuario. Also uote that the amplitude of (6.1) is 
mostly iudepeudeut of whether the iuhatiou model is of the large or small held type. 


• As the mass of the clock held becomes much larger thau the Hubble parameter if, all 
the amplitudes of the clock siguals above are expected to be suppressed by a Boltzmauu 
factor. This factor should be roughly of the form ~ exp(—7rm/if) as iu |^, aud the 
details are uuder iuvestigatiou. Therefore, for the observatioual purpose, the most 
importaut staudard clock helds should be those with mass larger thau 3if/2 but still 
of order if. luterestiugly, scalar helds with mass of order if are geuerically expected 
iu iuhatiou models iu supergravity. 


As we cau see, the amplitude of the clock siguals spaus a wide rauge aud is very model- 
depeudeut. While iu mauy cases the siguals are small, some of them, iucludiug some uatural 
oues iu the supergravity model-buildiug, may fall iuto the observable rauges of the CMB, 
LSS, aud 21 cm experimeuts. It would be iuterestiug to study the observability of the 
siguals iu future experimeuts. (See Ref. 19,20 for a forecast ou the uou-oscillatory siguals 
iu quasi-siugle held iuhatiou models.) We also uote that, due to the special shape-depeudeut 
oscillatiou patterus with variables such as frequeucy aud phase, these siguals are most likely 
orthogoual to the existiug popular templates used iu uou-Gaussiauity data aualysis. So 
dedicated templates aud data aualyses are uecessary to coustraiu or discover them, for 


example, by applyiug the method of modal decompositiou 33 


7 Conclusions and discussions 


Gosmological observables that cau be used to model-iudepeudeutly distiuguish the iuhatiou 
sceuario from possible alteruative sceuarios are very rare aud au importaut subject of theo¬ 
retical aud experimeutal study. Iu this paper, motivated by the developmeut iu the classical 
primordial staudard clock models (fll^ aud quasi-siugle held iuhatiou models [TS, 18 , we 
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pointed out and studied a novel type of such signals generically existing in the primordial 
universe models, independent of the well-known tensor mode approach. These signals are 
generated by quantum fluctuating massive helds, and are the quantum mechanical realiza¬ 
tion of the idea of the primordial standard clocks. The ticks of the massive held oscillations 
get imprinted as oscillatory features in the shape-dependence of various primordial non- 
Gaussianities, which directly encode the time-dependence of the scale factor of the primor¬ 
dial universe a{t). We have demonstrated how this is achieved through qualitative arguments 
as well as detailed examples in both inhation and alternative-to-inhation scenarios. In the 
inhation example, we are able to derive several exact results. In the alternative-to-inhation 
examples, we are able to use the universal features of the standard clocks to derive some 
general results. 

Here to highlight some key properties, we summarize the main differences between the 
classical standard clocks, studied in previous works, and the quantum standard clocks of this 
work. 


• Classical standard clocks exist only if some massive helds get excited classically by cer¬ 
tain background sharp features. In contrast, the quantum standard clocks generically 
exist in realistic primordial universe models because massive helds always huctuate 
quantum mechanically in time-dependent background. Depending on the structures of 
couplings in detailed models, the signals may show up in diherent types of three-point 
and/or higher order correlation functions, with the universal gravitational coupling 
case being the minimum case. However, the amplitude of the signals is highly model- 
dependent. This situation is much like that for the tensor mode. 

• Due to the generation mechanism, the simplest classical standard clock models predict 
two diherent kinds of signals that are smoothly connected to each other, namely the 
sharp feature signal and clock signal] while the simplest quantum standard clock models 
only have the clock signal. 

• Classical standard clocks break the scale invariance. The clock signals manifest as 
scale-dependent oscillatory features in various correlation functions including the power 
spectrum and non-Gaussianities. In contrast, the quantum standard clock itself is 
not a source of scale-invariance breaking p] because it appears persistently due to its 
quantum nature. So, the clock signals are absent in the power spectrum, and are only 
present in non-Gaussianities and manifest as shape-dependent oscillatory features. 


There are many important directions that need to be studied for the quantum primordial 
standard clocks. For example. 


This does not include the possible scale-dependence which may be introduced by the primordial universe 
model itself when clock fields are ignored, especially in the alternative-to-inhation scenarios. Relatedly we 


note that, for alternative-to-inhation scenarios, if we rescale all momenta in (3.15) and (3.16), the phases 
may not be invariant at the leading order due to the possible scale-dependence in Hk^ and Such a 

dependence does not encode the clock signal. 
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• We have assumed that the mass of the clock held remains constant. Time-dependent 
mass terms are naturally induced in the background where there is no shift symmetry 
in time. If the mass of the held is much larger than the horizon scale, this time- 
dependent contribution to the mass is negligible. This is usually the case for the 
classical standard clocks, in which the amplitude of the clock signal increases as the 
mass of clock increases [^[^. However, for the quantum standard clock models, the 
signal is exponentially suppressed by a Boltzmann factor for heavy clock helds, so the 
observable range is where the mass of the clock held is not too much larger than the 
horizon scale. The time-dependence in the mass term mentioned above can become 
important in this case, and it is interesting to study how it ahects the predictions. 


We have focused on the shape of non-Gaussianities, but not their amplitudes. The 
main point is to show the general existence and physical signihcance of the quantum 
standard clock. Nonetheless, the amplitude is critical for the observability of such 
signals, and as discussed it is highly model-dependent. A promising example is the 
example of quasi-single held inhation model studied in 13 . It will be interesting to 


extend the results in 13 to m > 3H/2 and compute explicitly the clock signal in this 
model. This could provide a concrete inhation model example in which the quantum 
standard clock signal can be observably large. 


We only considered one type of the three-point function with specihc interaction vertex 
depicted in Fig. The clock signals are also encoded in shapes of other types of three- 
point functions (with diherent vertices) and higher-order correlation functions, in terms 
of similar squeezed conhgurations or other special momentum conhgurations. This is 
another important direction. Related studies in quasi-single held inhation models can 
be found in 


15,18 


• In the analysis of the alternative-to-inhation scenarios, we have focused on generic 
properties of the quantum standard clocks that are insensitive to model-building de¬ 
tails. It would be interesting to study more explicit and concrete models. 


• In this work, we have emphasized the key theoretical predictions and physical signih¬ 
cance of the quantum standard clocks. An important future step is to survey diherent 
models, come up with simple templates for such signals and implement them in the 
actual data analyses. 
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A Clock signal and (non-)resonant integral 


As we explained in the main text, the ticks of the massive clock held oscillation get imprinted 
in the density perturbations and record the time-evolution of the background scale factor. 
This process is mostly done through the resonance mechanism [9 12 . In order for the 
resonance mechanism to take place, if the massive held contribution to the integrand behaves 
as ~ then the curvature scalar held contribution has to behave as ~ However, as 

we have seen in Sec. 4T, in the inhation case, the clock signal appears even if the two modes 
in certain integral do not have the required arrangement above and hence do not resonate 
with each other. In such a case, the resulting clock signal is suppressed by an exponentially 
small factor; but nonetheless it may still be important because such a non-resonant integral 
may carry an exponentially larger weight than the resonant one, due to some natural initial 
conditions. In Sec. A.l we illustrate this point in terms of simple integrals. 

We also study the same issue for the alternative-to-inhation models with 0 < p < 1. In 
Sec. A.2 through several representative examples we observe that, unlike the inhation case, 
resonance is required to generate the clock signal. However we do not have a general proof 
for arbitrary values of p within 0 < p < 1, and so cannot exclude the possible existence of 
special cases. This conclusion will be used in Sec. 


A.l Exponential inflation example 


The integral that we refer to in the inflation case is the term in the 2nd bracket of (|3.8|): 




(A.l) 


This definite integral is done exactly in Sec. 4T For example, in the very squeezed limit, 
^ 3/^1 ^ large mass limit, p S> 1, we have 




'Y) nit 


(A.2) 


We have chosen the massive held Vk to be originated from the BD vacuum. In the classical 
regime, this choice (after complex conjugation) leads to an oscillatory behavior that is a 
linear combination of the negative and positive frequency, {—x)~^^ and (—x)*^; but with 
very diherent weights - the former is a factor of larger than the latter. The positive 

. k-i + k2 

frequency component ~ (—x)®^ resonates with the massless mode ~ e* ''3 ^ in the integrand. 


giving a clock signal which is the flrst term in (A.2). Interestingly, although the negative 
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frequency component does not resonate, the convoluted integration with the massless mode 


contributes a term that also contains a similar clock signal (the 2nd term in (A.2)); there 


is much cancellation during this integration due to the absence of resonance, but this is 
compensated by the large weight of this component. In the end, the two clock signals in 


(A.2) are equally important, causing some phase shift to each other. We also note that, in 


some integrals such as Xnankei in (4.13), the massive held mode function differs from that in 


(A.l) by a complex conjugate. Then the component that resonates with the curvature scalar 


is different, and the non-resonant clock signal becomes negligible. 

To best illustrate the point, we look at the following toy integral. 


i±(p) = 


dre 


■^imt ^—iKr 


(A.3) 


where corresponds to the resonant integral, and X_ the non-resonant integral. This 


integral captures the essential behaviors of both oscillatory terms in (A.l), and for the 


purpose of this appendix, we ignore the other slowly varying functions in the integrand. 

In the case of exponential inhation that we consider in this paper, p —)■ cx), so the integral 
is in the following form 


X±(cx)) = 




(A.4) 


After using the ie prescription to make the integrals UV convergent, one can obtain the 
following exact results 


X+(cx>) = r(l - i^) (A.5) 

X_(oo) = +1). (A.6) 


It is clear that both resonant and non-resonant integrals give the clock signald^ However, 
the non-resonant integral is suppressed in comparison with the resonant one. This is indeed 
expected as the resonance enhances the amplitude, while the non-resonant integral is sup¬ 
pressed due to the highly oscillatory factors. To see this fact more clearly, it is useful to 
expand the exact results in large mass limit, i.e. p S> 1. One has 


X+(oo) + (A.7) 

X_(oo) (A.8) 

The extra factor of e~'^^ in X_ is responsible for the suppression of the non-resonance integral. 

^^For cases where p > 1 but still at the order of one, the situation is a little tricky. For example, we have 
checked that for p = 2 the non-resonant integral does not produce the clock signal while for p = 4 it does 
but with a different frequency from the resonant case. The more non-trivial case is p = 3 where the resonant 
and non-resonant integrals give exactly the same results and there are two oscillatory factors with different 
frequencies (though the scale dependence of the frequencies is the same, consistent with the clock signal). 
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A.2 Alternative to inflation (0 < p < 1) 


In this section, to study the same issue in the contracting scenarios (0 < p < 1), we use the 
same toy integral (A.3). Following the same convention in Sec. for 0 < p < 1 this integral 
becomes 


rO 1 

^±{p)= dr 


(A.9) 


We are unable to analytically integrate (A.9) with general p, so, in the following we compute 


this integral with representative p values for two types of contraction scenarios, respectively. 
The conclusion applies to these examples at least, although we do not have a general proof. 


A.2.1 matter contraction (p = 2/3) 
In this case we have 
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and 


I-(2/3) = 
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(A.ll) 


where 1 A 4 is the Generalized Hypergeometric function and Ai and Bi are Airy functions. 
Note that one can reproduce X_ from X+ by the analytic continuation u —)■ i>e~™. In order 
to see whether or not the results show the clock signal it is useful to expand the results for 
large arguments, where we expect characteristic oscillatory behavior in the clock signal: 


/ _ 2/y3/2 _ 2/y3/2 \ 

X+(2/3) ^ > oc (a/4 + 3i sin( ) + a/ 4 - 3i cos( )J (A. 12 ) 

(A.13) 


X_(2/3) 




— f> Zy/u _L - 

2(Az/)i/4 K ’ 


where we kept the leading order terms in the real and imaginary parts. In (A.12) we could not 


determine the precise prefactor directly from (A.IO). The above expansions show clearly that 


the resonant integral gives the clock signal expected from the saddle point approximation, 
while the non-resonant integral is just a decaying function in A. Note that the above 
approximate results are not related by the analytic continuation u —>■ That is, the 

expansion for large A and the analytic continuation do not commute. 
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A.2.2 A slowly contracting scenario (p=l/5) 


We next study the same integral in a slowly contracting scenario with p = 1/5. Again the 
integrals can be computed exactly and one has 


X+(l/5) 


5^/^ 9 / 9 7 19_ 2 3 4_ iK ^\ 
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(A.14) 


X_(l/5) 
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(A.15) 


Expanding the above exact results in large K limit one has 
X^(l/5) 2(1 gif/ ^ 

™ . 4(-l)5/353/4^Vr(l) 1 

77 + r(-A)r(A)r(i;)r(i) 7^9/4 ’ 


(A. 16 ) 

(A.17) 


where we have kept the leading order terms in the imaginary and real parts. Again, only the 
resonance integral is oscillating and the frequency of oscillation is consistent with the clock 
signal obtained by using the saddle point approximation. 

We have checked that the same conclusion applies to the p = 1/10 case. 
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